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Abstract 

Both particle physics and the 1890s Seeliger-Neumann modifica- 
tion of Newtonian gravity suggest considering a "mass term" for grav- 
ity, yielding a finite range due to an exponentially decaying Yukawa 
potential. Unlike Nordstrom's "massless" theory, massive scalar gravi- 
ties are strictly Special Relativistic, being invariant under the Poincare 
group but not the conformal group. Geometry is a poor guide to un- 
derstanding massive scalar gravities: matter sees a conformally flat 
metric, but gravity also sees the rest of the flat metric, barely, in the 
mass term. Infinitely many theories exhibit this bimetric 'geometry,' 
all with the total stress-energy's trace as source. All are new except 
the Freund-Nambu theory. The smooth massless limit indicates un- 
derdetermination of theories by data between massless and massive 
scalar gravities. The ease of accommodating electrons, protons and 
other fermions using density-weighted Ogievetsky-Polubarinov spinors 
in scalar gravity is noted. 

keywords: scalar gravity, Klein-Gordon equation, Nordstrom, massive, 
spinor, scalar density 
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1 Introduction 



In the 1910s Nordstrom proposed a theory of gravity that met the stric- 
tures of Special Relativity [H [21 E] by having, at least, Lorentz transfor- 
mations as well as space- and time-translations as symmetries, as well as 
displaying retarded action through a field medium, as opposed to Newto- 
nian instantaneous action at a distance. Nordstrom's scalar gravity was a 
serious competitor to Einstein's program for some years during the middle 
1910s. Neglecting time dependence and nonlinearity, it gives Poisson's equa- 
tion just as Newton's theory does. Nordstrom's theory was eclipsed first by 
the theoretical brilliance of Einstein's much more daring project, and then 
by the empirical success of Einstein's theory in the bending of light, a result 
inconsistent with Nordstrom's theory. While representing gravity by a scalar 
field is no longer a viable physical proposal, it is interesting to fill a hole left 
by the abandonment of Nordstrom's scalar gravitational theory caused by 
Einstein's inventing General Relativity (GR) so soon. 

Developments in group theory as applied to quantum mechanics, such 
as by Wigner jl], classified all possible fields in terms of the Lorentz group 
with various masses and various spins. In the late 1930s Pauli and Fierz 
found that the theory of a non-interacting massless spin 2 (symmetric tensor) 
field in Minkowski space-time was just the linear approximation of Einstein's 
GR El U\. Tonnelat and some author authors associated with de Broglie 
pursued massive spin 2 theories [5]. Nordstrom's theory of a long-range scalar 
field is, in this particle physics terminology, a theory of a massless spin field. 

Both the precedent of particle physics in the 1930s and the work by 
Seeliger and Neumann in the 1890s in giving Newtonian gravity a finite 
range [HI QUI EH E21 ESI EH ES] show the appropriateness of considering the 
possibility of a finite range for the gravitational potential. A finite range 
corresponds in field theory to a mass term, a term in the field equation 
that is linear and algebraic in the potential; the corresponding term in the 
Lagrangian density, a sort of potential energy, is quadratic. 

Einstein briefly entertained in his 1917 paper on the cosmological constant 
what is in effect a massive scalar gravitational theory 

as a foil for what is to follow. In place of Poisson's equation we 
write 

V 2 0-A0 = 4vtkp ... (2) 

where A denotes a universal constant. [HI p. 179] 
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Thus Einstein in effect contemplated a theory of the sort that, in light of 
later quantum mechanics terminology, one might call a theory of gravity 
using a massive scalar field, with A equaling the square of the scalar graviton 
mass in relativistic units with Planck's constant and the speed of light set 
to 1. Relativistic massive scalar fields in the absence of interacting satisfy 
the Klein-Gordon equation, but interpreting the field as gravity introduces 
interactions, including self-interaction and hence nonlinearity. 

However, Einstein promptly drew a widely followed analogy to his cosmo- 
logical constant — I suppress references to spare the guilty — but this analogy 
is erroneous [13 [18l [T9l [20l [2U [221 HH I2S1 H5]. Thus Einstein obscured for 
himself and others the deep conceptual issues raised by the mass term. The 
cosmological constant, having a zeroth order term in the field equations, is 
analogous to the scalar equation 

V 2 - A(l + 4>) = 4-KKp; 

the strange term —A • 1 will tend to dominate over the intended — \<p term. 

Two papers from c. 1970 provide a partial exception to the remarkable 
silence about the possibility of relativistic massive scalar gravity theories. 
One by Freund and Nambu [23] writes down equations that formally could 
be read as applying to massive scalar gravity, but they do not consider that 
application. Deser and Halpern [25] soon called attention to the identity of 
the Freund-Nambu field equations (less the mass term, though that omission 
is not mentioned!) with Nordstrom's theory. While the alert reader could 
notice that Freund and Nambu had in effect provided field equations for a 
massive variant of Nordstrom's theory, apparently no one has ever managed 
to comment on that fact, still less to discuss its significance. A more recent 
paper by H. Dehnen and R. Frommert reinvented massive scalar gravity 
in a similar way [2S], but with an unfortunate and apparently unnoticed 
restriction on the allowed matter field content such that standard scalar 
fields are inadmissible. 

Massive scalar gravities, if the mass is sufficiently small, fit the data as 
well as does Nordstrom's theory, as a consequence of the smoothness of the 
limit of a massive scalar field theory as the mass goes to zero [23 P- 246]. 
Thus there is a problem of underdetermination between the massless theory 
and its massive variants for sufficiently small masses [28J. The analog of this 
instance of underdetermination was already clearly understood by Seeliger 
in the 1890s. He wrote (as translated by John Norton) that Newton's law 
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was "a purely empirical formula and assuming its exactness would be a new 
hypothesis supported by nothing." [2S1 E] While that claim might be a 
bit strong, in that Newton's law had virtues that not every rival formula 
empirically viable in the 1890s had, a certain kind of exponentially decaying 
formula of Neumann and Seeliger was also associated with an appropriate 
differential equation [TU]. 

It is well known that Nordstrom's theory does not bend light [SU]- That 
is an immediate consequence of the conformal flatness of the metric in Nord- 
strom's theory in geometrical form [3T| [32] and the conformal invariance of 
Maxwell's electromagnetism [33]: space-time is flat in Nordstrom's theory 
except for the volume element, but light doesn't see the volume element in 
Maxwell's theory. While scalar gravity is a museum piece as far as theoret- 
ical physics is concerned — at least as far as the dominant gravitational field 
as concerned — it remains a useful test bed theory for analogous phenomena 
for which the details in General Relativity are much more complicated tech- 
nically or might play a secondary role in gravitation theory [3U |35l [361 137] . 
Scalar gravity sheds light on some fundamental issues in space-time theory 
as well, not least by setting a precedent that could be entertained for massive 
tensor theories. 

2 Massive Scalar Gravities: Relatives of 
Nordstrom's Theory 

Here I shall give a suitable Lagrangian density in a form adapted to the 
derivation of the massive variants to be introduced shortly. The Einstein- 
Fokker geometrization [2D E2] suggests a useful set of variables to use. One 
can isolate the conformal structure (the null cones) out of a metric by taking 
the part with determinant of —1; for a flat metric T]^ one can call the resulting 
conformal metric density fj^, a tensor density of weight — | in four space- 
time dimensions [38]. Let fj be a (positive) scalar density of arbitrary nonzero 
weight w, related to r/^ by ^/—det(rj^) = f\™. (The expression \J —det^^) 
is often written as y/—rj-) Thus fj (or its wth root) governs volumes, at least 
volumes that are not distorted by gravity. Let the gravitational potential be 
represented by a potential 7, also of density weight w. Then one can define a 
new effective volume element by g =def fj + kj. Thus far it is unclear whether 
g actually determines the volume of anything, but the derivation shows that 
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it determines all volumes for the matter dynamics and all (in the massless 
case) or most (in the massive cases) for the gravitational dynamics. It turns 
out that k 2 = 64irGw 2 ] the sign of k does not matter much, but will be 
chosen to match that of w to maximize continuity. Neglecting terms that 
do not contribute to the field equations or that disappear with the choice 
of appropriate coordinates, one can take the purely gravitational part of the 
Lagrangian density as 



here is the coordinate derivative, while rf u is the inverse of f)^ u . This re- 
sult will be derived below. It is essential to use scalar densities in order to 
obtain the trace of the stress-energy tensor so readily, which then permits 
combining the gravitational potential with the background volume element 
by an additive field redefinition g =def fj + kj. The basic postulate is univer- 
sal coupling, that the full field equations are obtained by taking the free field 
equations and adding in the trace of the total stress-energy tensor (including 
gravitational energy-momentum). The universal coupling principle (initially 
using a brute-force direct construction of the stress-energy tensor) was em- 
ployed by Einstein in his supposedly unsuccessful Entwurf physical strategy 
for finding his field equations [39], and was later brought to successful com- 
pletion in finding Einstein's equations using higher mathematical technology 
[30] . One can show that the trace of the stress-energy tensor is given by 
taking the Euler-Lagrange derivative of the Lagrangian density with respect 
to the volume element (perhaps raised to some power), such as 



The same theory results from all values of w. Cf ree is a rather standard 
quadratic expression yielding the Klein-Gordon equation. The material part 
of the Lagrangian density can be written in terms of matter fields u and 
the conformally flat metric built from fj^ u and g; fj does not appear on its 
own. Both the gravitational and the material parts of the Lagrangian density 
therefore depend only on a conformally flat metric. 
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While Nordstrom's theory has been derived in terms of universal coupling 
previously J3UJ I2U 125] , the approach outlined here admits a ready general- 
ization yielding a distinct massive scalar gravity for every value of w, only 
one of which has been found before. The kinetic and matter terms are of 
course just those of Nordstrom's theory Given the nonlinearity of the the- 
ories, there are different choices of field variables, often nonlinearly related, 
that are especially convenient for one purpose or another; thus comparison re- 
quires choosing some common set of fields. It is convenient to use \f—rj = fj™ 
and y/—g = g™\ here i] and g (without the ~) are the determinants of the 
metrics rj^ and g^ v as usual. In the mass term only, one has sJ—T) appearing 
in its own right, the fact around which much of the philosophical interest of 
the theories. For any real w (including w — 1 and w = by l'Hopital's rule), 
a universally coupled massive variant of Nordstrom's theory is given by 
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One can express this mass term as a quadratic term in the potential and, 
typically, a series of higher powers using the expansion 
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-rj W + 8wV kG^y, 



where 7 is the gravitational potential; the case w = requires taking the — th 
root of this equation and taking the limit, giving an exponential function. In 
all cases the result is 
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This one-parameter family of theories closely resembles the 2-parameter fam- 
ily Ogievetsky-Polubarinov family of massive tensor theories jlOj HI] , which 
can also be derived in a similar fashion jl2]. The case w = |, which conve- 
niently terminates at quadratic order, is the Freund-Nambu theory [21]. It 
is also useful to make a series expansion for all the theories using the special 
w = exponential variables, as will appear below. 
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3 Massive Scalar Gravities Are Strictly Spe- 
cial Relativistic 



Features of Nordstrom' s scalar gravity are said to have shown that even the 
simplest and most conservative relativistic field theory of gravitation burst 
the bounds of Special Relativity (SR) [HI p. 179] [H P- 414 ]- Relativistic 
gravity couldn't be merely special relativistic, according to this claim. Nord- 
strom's theory indeed has a merely conformally flat space-time geometry [31] , 
which one can write as 



where rj^ (with determinant —1) determines the light cones just as if for a flat 
metric in SR. Thus Nordstrom's theory is invariant under the 15-parameter 
conformal group, a larger group than the usual 10-parameter Poincare group 
of Special Relativity. 

By contrast the massive variants of Nordstrom's theory are invariant only 
under the 10-parameter Poincare group standard in SR and thus are special 
relativistic in the strict sense. The mass term breaks the conformal symmetry. 
It is therefore false that relativistic gravitation could not have fit within the 
confines of Special Relativity. While it is true that no phenomena required 
the mass term, it was possible that the mere smallness of the mass parameter 
explained its empirical obscurity, as Seeliger had already proposed in the 
Newtonian case. 



4 Derivation of (Massless) Nordstrom Scalar 
Gravity 



Having summarized the results above, I turn to their derivation. While 
Nordstrom's theory has been derived in terms of universal coupling previously 
[30~t Ell EE] , the approach outlined here admits a ready generalization yielding 
a distinct massive scalar gravity for every value of w. All of these theories are 
new, with the exception of the Freund-Nambu theory. The use of a scalar 
density as the field variable permits a convenient additive change of variables 
and simple coupling to the trace of the total stress-energy tensor. A scalar 
density under arbitrary coordinate transformations is of course a scalar under 
Lorentz transformations (with two options regarding transformations with 
negative determinant [35], corresponding to the scalar and pseudoscalar of 
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particle physics). It will turn out that different density weights conveniently 
give different massive scalar gravities, much as in the tensor case [101 H2] • 

4.1 Tensor Densities and Irreducible Parts of the Met- 
ric and Stress Tensor 

In addition to the background geometrical variables and the gravitational 
potential, there are also other matter fields, which I denote collectively by u. 
The derivation that follows makes use of the Rosenfeld-type metrical defini- 
tion of the stress-energy tensor [15J EDI El SSI US]- In such an approach the 
flatness of the metric is momentarily relaxed while a functional derivative 
with respect to it is taken; then flatness is restored. It will be most help- 
ful to use not the metric tensor rj^, its inverse, or some densitized relative 
thereof, but rather two irreducible geometrical objects that together build up 
the metric tensor. (The use of irreducible geometric objects is also impor- 
tant in analyzing the Anderson- Friedman geometric objects program [5U].) 
The conformal metric tensor density fj^ (with determinant of —1) of weight 
— \ determines the null cone structure, which is untouched by gravitation in 
scalar gravity theories. The remainder of the flat metric tensor is supplied 
by a scalar density fj of nonzero weight w, which quantity we may take to be 
positive^ The flat metric tensor is built up as rj^ = Tj^uiv)^ ■ Note that w 
can be any real nonzero number; the use of a single tilde above the variable 
name reminds us that the field is a density, but gives no hint about which 
weight it has. The usual torsion-free covariant derivative compatible with 
r]^ v is written as d a . The metric-compatibility condition d a r]^ u = yields 
d a fj = and dafj^ = as well. One should be careful to use the correct 
form of the covariant derivative for tensor densities, so it is worth recalling 
the forms of their covariant and Lie derivatives. A (1,1) density 0? of weight 
w is representative. The Lie derivative is given by [51j [38j E2] 

£^1 = er^ -0?e% , (7) 

where the , \x denotes partial differentiation with respect to local coordinates 
x^. The ^-covariant derivative is given by [5TJ [381 E2] 

d,^ = <^ - cr^ - w^v^. (8) 

1 Choosing fj > in all coordinate systems indicates which of the various subtypes of 
density [45] is intended. 
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Here Tp are the Christoffel symbols for rj^. Once the curved metric 
is defined, the analogous g-covariant derivative V with Christoffel symbols 
follows. For scalar densities, the relevant piece is the new term with 
the coefficient ±w. The formulas for Lie and covariant differentiation follow 
[53| 151] from the coordinate transformation law for scalar densities: under a 
change of local coordinates from an unprimed set x^ to a primed set x u , the 
density's component behaves as 



(Ox 
to 



(9) 



The primes perhaps are opposite where one might have expected, but this 
is the usual convention [381 113 EE] , though some authors, especially but not 
only Russian, define weight in the opposite fashion. For the other kind of 
scalar densities, which can change sign under coordinate transformations [45J, 
the same Lie and covariant derivative formulas follow, because only behavior 
for infinitesimal transformations is relevant. 

For any action S invariant under arbitrary coordinate transformations, 
one can derive a metrical stress-energy tensor. It is convenient to break the 
flat metric up into its irreducible parts, the conformal metric density f]^ that 
fixes the null cones and the volume-related scalar density fj [19]. One can 
show that 
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Making an infinitesimal coordinate transformation described by the vector 
field gives 
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where BT is some boundary term of no interest for present purposes. Be- 
cause S is a scalar, 5S = 0. Integration by parts pulls all the derivatives off 
£ M at the cost of more boundary terms. Choosing £ M to vanish at the bound- 
ary annihilates all the boundary terms, while using its arbitrariness removes 
the integration, leaving the integrand to vanish. Going 'on-shell' by using 
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The stress-energy tensor here is broken up into traceless and trace parts. 
For a scalar theory, the latter will represent the source for gravity. This 
stress-energy tensor contains contributions from both matter u and gravity 
7. If Maxwell's electromagnetism is included among the matter fields, then 
it couples only to 77^, not 77: rather than introducing conformal rescalings 
and showing that nothing interesting changes [331 ES], this use of densities is 
the most direct way to see the theory's conformal invariance, from which the 
theory's failure to bend light follows immediately: electromagnetic radiation 
in the absence of charges doesn't know that it isn't in flat space-time because 
the volume element disappears entirely. 

4.2 Spinors 

Foundational questions about space-time not infrequently are addressed as 
though there were no such thing as protons, electrons, neutrons, or other 
fermions, which are represented classically by spinor fields, or as if spinor 
fields introduced no additional issues. But consider the supposedly triv- 
ial possibility of representing special relativity using arbitrary coordinates. 
While mere tensor calculus is needed for bosonic fields, the issue is more 
complicated for spinor fields. In the interest of overcoming this unjustifiable 
neglect of spinor fields, I point out that spinorial matter can be included 
effortlessly in the universally coupled massive gravities considered here. On 
account of breaking up the metric and stress tensor into their irreducible 
pieces, including spinors requires no work at all, because the spinor (if mass- 
less) does not notice scalar gravity. To see this effortlessly and yet rigorously, 
one can use the Ogievetsky-Polubarinov spinor formalism [lOl [561 E3 158] - 
thereby avoiding a tetrad in favor of the metric (which is possible by con- 
struction, contrary to widely held belief); the formalism is a bit like the tetrad 
formalism in the symmetric gauge, but is conceptually independent. The con- 
formal covariance of the massless Dirac equation [551 f59j is well known. But 
what is rarely if ever noticed is that one must and can use density-weighted 
spinors to achieve conformal invariance, with the volume element dropping 
out altogether — much as in Maxwell's electromagnetism, but the details are 
more difficult and less familiar and involve derivatives of the conformal part 
of the metric [60] . The trace of the stress-energy tensor comes from || , but 
fj is simply absent from the suitably weighted spinor's kinetic term. That ex- 
pression depends, in a highly nonlinear way, only on the components of fj^, 
which determines the null cones: 9 components, not the 10 of the metric, 15 
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of a unimodular conformal tetrad, or 16 of a tetrad. The appropriate spinor 
has density weight | in 4 space-time dimensions or in n space-time di- 
mensions [60] . The spinor, if massless, does not notice scalar gravity at all; 
the metrical stress-energy tensor has vanishing trace even off-shell. More fa- 
miliar routes to this conclusion of vanishing trace of the metric stress tensor 
are less direct [BT| W2\ IB15] . The Ogievetsky-Polubarinov treatment of spinors 
has also avoided a spurious counterexample to the Anderson-Friedman ab- 
solute objects program [50] for understanding the difference between merely 
formal general covariance and the substantive kind that is supposed to be a 
novel feature of General Relativity. 

4.3 Universal Coupling 

Let us assume that the free field action S/[7, u, fj^, fj] (with vanishing New- 
ton's constant G) is known; it is given below, apart from the unspecified 
matter fields. The full action 5* should reduce to Sf in the limit of vanishing 
G. The task at hand is to derive the full action S for the theory with nonzero 
gravitational interaction. Lorentz covariance requires a source generalizing 
the mass density in a Lorentz-invariant way, so the trace of the stress-energy 
tensor is a natural choice. This choice is perhaps not compulsory, unlike 
the tensor case where free field gauge invariance necessitates that any source 
used by a divergenceless symmetric rank 2 tensor, of which there is only one 
physically significant example at hand. But it is a very natural choice. 

Letting S be specialized once more to the interacting scalar gravitational 
theory that we seek, we can postulate that the gravitational free field equation 
is modified by the introduction of a source term that is basically the trace of 
the stress-energy tensor: 

5S _ SSf +k ^S 
#7 8q 5fj ' 

where k is a coupling constant related to G and perhaps the density weight 
w in ways that will be ascertained later. In anticipation of a change to 
bimetric variables, the |7 notation has been added to emphasize that the 
other independent variable here besides fj is 7. 

The change to bimetric variables^ involves the definition 

2 For the tensor case, the new variables involve two metric tensors or the like 30, 49, 42 . 
For the scalar case, only the scalar density portions undergo any redefinition, because the 
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Equating coefficients for variations of fj and 7 relates the functional deriva- 
tives as follows: 

5S t „ 5S 6S.„ 
of) og of) 

— -k— 

$1 $g" 

The second equation shows that the new field g has an Euler-Lagrange equa- 
tion equivalent to that of the potential 7. The first equation shows that the 
trace of the stress-energy tensor splits into one piece that vanishes on-shell 
and one that does not. Using these equations in the postulated equation for 
universal coupling gives 

5Sf , 8S 
0=-f+k—g. 
07 or) 

So far little has been said about the detailed form of Sf. The most natural 
choice is given by the Lagrangian density 

apart from a mass term for 7 which will be included later (and a free field 
term for matter u which does not contain 7 and so does not contribute to 
the derivation). This choice gives the usual wave equation. The metrical 
signature — h ++ is employed. 

To satisfy the universal coupling identity in bimetric guise, it is convenient 
[30l [38j H9] to split the full (but unknown) action S into a piece Si[g, u, f)^ u ] 

(without explicit dependence on fj) and another piece S 2 that, perhaps among 

ss 

other things, cancels the term -jJ-. For Si it is natural to build a conformally 
flat metric 

9vu = f)^(g)^. 
Then it is natural to choose the Hilbert-like expression 

Si = c J d 4 Xy/gR[g] + S maU er[gfMu,u\. 

(A cosmological constant term J d^x^J—g is also available if desired.) One 
can show that 

S 2 = ^ J d 4 xR[r)]f)- l+ -l + J d 4 xd^ 



gravitational field lias no tensor piece to combine with . Nonetheless the term "bimetric 
variables" is handy. 
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does the job of accommodating Sf. The first piece Si involving the Ricci 
scalar for r)^ does the work here. The second piece S 2 is simply a bound- 
ary term, which is deposited into S 2 rather than elsewhere for convenience. 
The boundary term can be chosen to remove the second derivatives from the 
Hilbert-like term y/gR\g\ in S±. One can also include a pure volume term 
J d 4 Xy/—rj. Note that S 2 contributes nothing to the field equations; its pur- 
pose is to contribute to the Rosenfeld metric stress-energy tensor only. (Re- 
call that flatness of the background is relaxed briefly in taking the functional 
derivative and then restored.) The total action for the massless case is thus a 
piece S\ describing an effective conformally flat geometry and a piece S 2 that 
does not affect the field equations. Universal coupling has completely clothed 
the background volume element with the gravitational potential, leaving only 
their sum as observable. This is an amusingly strict realization of Einstein's 
Poincare-inspired dictum that only the epistemological sum of gravity and 
physics is observable [61]. By requiring the usual normalization for the free 
gravitational field's kinetic term to lowest order, one infers that 

Aw 2 

c= ~w 

In comparison to the free gravitational Lagrangian density 

- \r u v^- 2 m)d u j, (10) 

the interacting theory has the corresponding expression (apart from terms 
not affecting the equations of motion) 

- \r u 9^- 2 (d^)d^ } (11) 

while the matter fields see the conformally flat effective metric rather than 
the flat background metric. Thus universal coupling of gravity to the trace 
of the total stress-energy tensor yield's Nordstrom's theory. In the massless 
case, the same theory obtains for every (nonzero) density weight w of the 
gravitational potential. The case w = without a graviton mass term was 
handled by Kraichnan [30J using an multiplicative exponential field redefini- 
tion, rather than an additive one, and yields Nordstrom's theory as well. 
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5 Derivation of Massive Gravities 



Making the gravitational potential 7 massive implies that the free gravita- 
tional potential obeys the Klein-Gordon equation, not the wave equation (or 
Laplace's equation in the static case). It is therefore necessary to add a term 
quadratic in the gravitational potential to the free gravitational field La- 
grangian density. Going through the universal coupling derivation, one then 
fields a crucial new term in which the flat metric's volume element remains 
essentially in the theory. One also crucially uses the cosmological constant 
term; the cosmological constant term and the new term with the flat volume 
element have equal and opposite terms linear in the gravitational potential in 
the total action, a crucial cancellation that removes the odd behavior (from 
a field-theoretic point of view |20j) of the cosmological constant. Then the 
pure volume term enters the action to cancel out the zeroth order parts of 
both the cosmological constant and the new term essentially involving the 
flat volume element; this last cancellation is largely a matter of good book- 
keeping, letting the action vanish for flat space-time but not affecting the 
field equations. Thus the interacting theory has a term quadratic in the 
gravitational potential (a mass term) and, except in one special case, one or 
more higher powers (possibly infinitely many) in the gravitational potential. 
While the derivation is carried out using different fields for different values of 
density weight w, making comparison mildly nontrivial, one can show (such 
as by using the two metrics only) that the theories obtained are all distinct. 
One can also show that the resulting one-parameter family of mass terms is 
related in the expected way to the two-parameter family of massive tensor 
gravities obtained some time ago by Ogievetsky and Polubarinov. (It turns 
out that the mass term in [40] contains a typographical error absent in the 
less well known summary |41j.) 

One expects that the mass term for a free field be quadratic in the po- 
tential and lack derivatives. The free field action Sf = S/o + S/ m is now 
assumed to have two parts: a (mostly kinetic) part Sjq that as in the mass- 
less case above, and an algebraic mass term Sf m that is quadratic. We seek a 
full universally coupled theory with an action 5* that has two corresponding 
parts. The two parts of S = So + S ms are the familiar part Sq (yielding 
the Einstein tensor, the matter action, a cosmological constant, and a ze- 
roth order 4-volume term) and the part S ms that essentially contains the 
background volume element fj and also has a zeroth order 4-volume term. 
As it turns out, the mass term is built out of both of the algebraic part of 
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So (the cosmological constant and 4-volume term) and the purely algebraic 
term S ms . In comparison to derivations using the canonical stress-energy ten- 
sor [2U |2Q], the metric stress-energy tensor is much cleaner, in some respects 
more illuminating, but in some respects less transparent. The canonical ten- 
sor derivation is noticeably simpler in the known w — ^ case for Freund and 
Nambu than in the other cases, for reasons that will be explained below. 
Again we postulate universal coupling in the form 

^-^l + k— (12) 
#7 <Py 5fj 

Changing to the bimetric variables g and fj implies, as before, that 

Now we introduce the relations Sf = Sfo + Sf m and S = So + S ms to treat 
separately the pieces that existed in the massless case from the innovations 
of the massive case. Thus 

SSm SSf m 5S .„ j8S ms .„ , . 

-jp- H — p— = —k-^r\g - k—zz-\g. 14) 
07 07 or) or) 

Given the assumption that the new terms S/ m and S ms correspond, this 
equation separates into the familiar part = — ^~fjH<7 as before and the 
new part 

ms 1 ~ 
^7 5fj 

So is given by So = Sifg^w] + S2 as in the massless case. Once again, we 
choose the simplest case and get the Hilbert action along with a cosmological 
constant, with matter coupled only to the curved metric, along with various 
terms that do not affect the equations of motion. Assuming the free field 
mass term to be quadratic in the gravitational potential, 




Cf m = — —1 yf^V , (15) 
its contribution to the field equation is 



57 

Changing to the bimetric variables gives 

m 2 ~~±-2 m 2 ~.L-i ,SS ms 
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5.1 Cases with w ^ 1, w ^ 

The goal for S ms is to obtain a mass term, so one can omit u and fj^ from 
the 'constant' of integration, leaving a function of sf—g only — which must 
be linear in order to be a scalar density of weight 1. Thus 

n f ,4 f a~±. wm ? i- 1 wm 2 „1_ 

S ms = Jdx (Ag^ + pp^N* - 

as long as w ^ 1. (The case w = 1 yields a theory as well, but must be 
treated separately.) Requiring S ms to vanish to zeroth order in 7 yields 

2 2 
w m 

A 



{w - l)fc 2 

Requiring to vanish to first order, which is important for the field equa- 
tions, gives nothing new. For the second and higher order terms, the binomial 
series expansion yields 

2 / 7 2 [l-2w}k^ 
*-ms = m y—V q h • 

For the free-field limit A; — ?• this expression reduces to the expected 
quadratic expression. For the special case w = \, which turns out to be the 
Freund-Nambu theory (the only case previously obtained), the mass term 
contains no interaction part; the free mass term, quadratic in a weight | 
potential, does not depend on the volume element in order to make a covari- 
ant (that is, weight 1) Lagrangian density and so contributes nothing to the 
metric stress energy tensor. (This perturbative expansion indicates nothing 
odd about the w = 1 case, though the above integration was not permissible 
in that case.) In terms of bimetric variables, the mass term for w 7^ 1 is 

w 2 m? wm 2 1 wm 2 _j_ 

[w — l)k z k z (l — w) k l 

The factor often can be treated using a binomial series expansion; the 
series converges for |/c7/r/| < 1 |65j. (A strong- field expansion is also possible, 
but will not be employed here.) One can show that the binomial series 
expansion for the theory labeled by w in terms of the weight w variables is 
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Here the expression 

-2)! 

is shorthand for (— — 2)(— — 3) • • • (— — j + 1); one need not make sense of 
the numerator and denominator separately in terms of Gamma functions, 
though one could do so. This form is clearly well behaved in the vicinity of 
w — 1, so one can find the limit as w — > 1 to be 

Note that this series expression of the theories leaves them not readily com- 
mensurable, due to the use of a different potential, bearing different relations 
to the more physically meaningful y/^g and y/—rj, for each value of w. The 
use a w-specific potential in the derivation of each theory is quite important 
in the context discovery, however. 

5.2 Case w = 1 

The case w — 1 can be considered now. The equation to integrate is 

m2 ~~-i m2 ,^s ms 

~T 9V -T = W lg - (19) 

Performing the integration introduces a logarithm: 

7T7.2 TH? 1 

Cms = -jj9[ ln (v) + f(g)} - -jjv- 

In the interest of getting a scalar density, one can set the 'constant' of inte- 
gration f(g) to be 

f(g) = a - ln(g) 

for some constant a. For w^l the quadratic and higher terms came from the 
formal cosmological constant term proportional to y/—g, but in this case that 
term is linear in the gravitational potential (in this choice of field variables), 
and hence serves to cancel the noxious linear part of the mass-yielding non- 
linear expression —^gln . This cancellation is a second service formed 
by the choice of —1 for the coefficient in /. (This remark should be taken 
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as merely heuristic, because use of nonlinear field redefinitions, such as re- 
expressing this mass term using a potential of a different density weight, or 
using the w = exponential field redefinition, would lead to a different sort 
of bookkeeping.) Requiring the zeroth order part of the mass term to vanish 
as well gives a — 1. Using the Taylor expansion ln(l + x) 
convergent for — 1 < x < 1 [66, p. 564], one obtains 



x 2 -r 3 



*-* r. 



m 



I 2 

2fj 



kj 3 
6fj 2 



matching the expansion above to this order. The full series expansion (for 
any w, with possible exception of w = — but that case will be vindicated 
shortly) can be shown to be 



C 



m 



-T] 



ms,w=l 



k 2 



E 

j=2 



-k- 
V 



(20) 



in agreement with the expression above for the limit of the 1 family in 
the limit as w — > 1. Thus the family of massive universally coupled scalar 
gravities is indeed continuous across w — 1, despite the need for special 
treatment of this case. 

One can also treat the case w = 1 using l'Hopital's rule for the indeter- 
minate form jj. One has 




64nG 



(21) 



where the formula for exponentials of a non-natural base introduces the log- 
arithms. This expression of course agrees with that given above. 
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5.3 Case w = 



The case w = is much more problematic, given the above bimetric field 
redefinition and the expression of universal coupling in terms of the derivative 
with respect to a volume element of some weight. The weight power of 
volume element is just 1, hardly a good field with respect to take a functional 
derivative. The additive field redefinition defining g appears to fail also. From 
the Newtonian limit it follows that 

k 2 = QAttGw 2 . 

To assess continuity of the field redefinition, one needs to know what happens 
to the meaning of 7 as w — > 0. By considering 

1 9wQ—w \J]w kw^fw)(J)—w w^f— w) ~ 1 ~\~ fjwk—w'J—w ~t~ k w ^f w fj— w 

near w — 0, one sees that the physical significance of the potential 7 does not 
jump discontinuously at w = as long as k and w have the same sign, which 
I choose to be positive for positive density weights w. Thus k = 8V irGw. 
But with k proportional to w, it appears that the bimetric field redefinition 

g w =fl + k(w)% 

(where the dependence on w has now been made explicit) reduces to 1 = 1 
for w — 0. The universal coupling postulate 

SS 5Sf , 8S 

— - = —^r- + k— (22) 
07 07 of) 

suffers not only from the meaninglessness of ||, but also from the vanishing 
of the coupling constant due to the linearity of k in w. 

While these problems seem rather disastrous, in fact they are all soluble. 
First we recall the series expansion above, with k now expressed in terms of 
w : 

C ms = -m ^(8 VnGY ^ j — , 

(23) 

which has well behaved and simple coefficients as w — > 0. It is natural to 
drop the tilde on 7 and set fj to 1 for w = 0, leaving the simple form 

OO A 
_ . , ryj 

,2 /—^X^(o./ZFi\3-2i 



i=2 J ' 
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This is clearly the sum of the quadratic and higher terms for the exponential 
function, so one infers 



C 



m 



ms,w=0 



-v 



64vrG 



-1 - 8VttG7 + exp{8V7rG~f) 



(25) 



It remains to find a meaningful and appropriate notion of universal coupling 
that permits, one hopes, the derivation of an expression equivalent to this 
w = series. 

The form of this series for w = suggests that one might think in terms 
of exponentials or logarithms to find a suitable field redefinition. While the 
linear redefinition 

g = fj + 8 V TiGw^i 

fails for w — 0, the ~th root 



v-g = g w 

makes sense for w 



(rj + 8v nGwy) ™ = 

also. The limit is 

\f—g = \/—T]exp(8V ttG^) 



7 



7? 1 + 8wV7rG-r 
r] 



(26) 



An exponential change of variables very much like this was already employed 
by Kraichnan, though without application to massive theories [30] . 

By writing the trace of the stress-energy tensor in two different ways, one 
can show that the problems of meaningless field variable and of vanishing 
coupling also can be resolved. The flat metric i]^ can be written as 



Thus one recalls that 



5S 1 5S 



Using this result in the postulate of universal coupling, the dependence on 
w cancels out, giving for w = 



— = -t 1 + 8wVnG— 

07 07 



6S 



2w Sr]^ u 



SS f A r— 5S 



5^/ 



6n 



(27) 
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which makes sense even for w — 0. It is convenient to choose the weight 1 
variable y/—rj, in terms of which universal coupling is 



SS 



SS f 8V^G SS 
1 + 



The exponential change of variables, while leaving fj^ u and u alone, gives 



SS 



7 



SS 



-9 



+ 



SS 



7] Sy^rj 
SS r- 



S~f S^-g 



-g- 



Installing these results in the universal coupling postulate yields 







SS 



#7 



S - + 8v / ^G v /3 ^ 



SS 



S^—T] 



\9, 



(28) 



(29) 



a result that is surprisingly indifferent to the non-additive form of the field 
redefinition. Letting the action S be a sum of Si + S2 from the massless case 
and S ms for the mass term, one has 



= - m 2 v /3^ 7 _|_ 8V7tGa 



SSr, 



-77- 



\9- 



(30) 



Making the change of variables in Sf as well yields 



= 



m 



-n 



In 



-77 



+ 8VttG, 



|0- 



Dividing by y/—r] and integrating gives 



(31) 



-647rG£ r 



[-v^^v^ + V^ln^r] - + /(#)] , 



where f(g) is a 'constant' of integration. To get a scalar action, the obvious 
choice is b^/—g for some constant b. Requiring the action to vanish to zeroth 
order yields b — 1; it vanishes to first order as well. The result is 



m 



64nG 



—\J—r\ln 



-9 



m y/—Tf 

64ttG 



-8VnG-f + exp(8ViiG-f) - 1 



(32) 
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which was already obtained above as the w — > limit of the series derived 
for the w 0. Thus the w = case in fact makes perfectly good sense and 
yields the theory that the w — > limit leads one to expect. If Kraichnan had 
considered massive scalar gravity, then he would have obtained the w = 
theory readily. 

One can also treat the case w = using l'Hopital's rule for the indeter- 
minate form jj. One has 



m 



lim 

™->o QAnG 



w — 1 



+ 



-1— w 



m 



lim 

™->o QAttG 



w(l — w) 



w \/—g — V~ 9 W y/~ V 1 w — (w — l)v^~^ 



w 



— w 



m 



lim 

w^o QAnG 



2w-l 



m 



64nG 



-rjlny/—rj + y/—rj), (33) 



in agreement with the formula given above. 

To sum up, while the w — 1 case needed some special treatment and the 
w = case needed a great deal of special treatment, every real value of w 
yields a (distinct) universally coupled massive scalar gravity. Thus we have 
found uncountably infinitely many massive scalar gravities, all derived by uni- 
versal coupling, that give finite-range rivals to Nordstrom's massless scalar 
theory. These theories all provide a relativistic embodiment of the Seeliger- 
Neumann finite-range modification of Newtonian gravity. There might be 
still other massive scalar gravities worthy of discovery, so no claim of ex- 
haustiveness is made. 



5.4 All Cases in w = Exponential Variables 

Above the infinite family of theories was presented both using the bimetric 
variables \f^g and \J—r\ and using a series expansion of each theory in its 
own adapted perturbative field j w . Recently it was found that the w = 
case suggests the relationship 

~g = exp(8V 7rG'-y)^—r); 
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this field 7 (with no~and no density weight) is a neutral, ecumenical choice for 
expressing all of the massive gravities in a commensurable fashion — as com- 
pared to the series expansions above, which use different fields for different 
theories. The result is 

m 2 y/=rj [we 8 ^ - e 8 ™ 7 ^ + 1 - w\ 
ms = ~64^G w(w - 1) 

for w ^ 0, 1; these special cases are readily handled by l'Hopital's rule. Us- 
ing the series expansion for the exponential function, which converges every- 
where, one has 

ms ' w ~ UtxG ^ j\ ^ 

j=2 J t=0 

I will not attempt to derive all the infinitely many theories using the w = 
exponential change of variables. While such a derivation must be possible in 
some sense, the premises might look contrived by virtue of the apparently 
non-linear form of some of the terms. Thus the role of the w-adapted field 
variables in the context of discovery is evident. They allow infinitely many 
derivations to succeed using a manifestly free field via a quadratic Lagrangian 
density, coupled to the total stress-energy tensor's trace, without powers of 
7 in the coefficients. 



6 Stability 

In the interest of avoiding runaway solutions due to a potential energy with no 
lower bound, one wants to investigate the behavior of the algebraic mass/self- 
interaction term. One might worry, for example, about theories such that this 
potential behaves like an odd polynomial (or worse) for large values (positive, 
or negative unless the singularity as y/—g — > matters-in fact it will prove 
helpful) of the gravitational field 7 or some relative thereof [HI]. While 
such strong values might invalidate the assumed validity of perturbation 
expansions sometimes made in this paper, a perturbative treatment at least 
suggests where trouble-spots might be found. For theories with the self- 
interaction potential behaving like an even polynomial (or certain kinds of 
infinite series), that sort of instability is not an issue, but correct physical 
interpretation requires checking whether 7 = or the like is the true vacuum 



23 



[57] . Checking these issues for all values of w would be a substantial task, 
but it is not difficult to check some interesting cases. Veneziano remarks 
that the Freund-Nambu theory is satisfactory on this count; I observe that 
for sufficiently negative values of the field there is a crushing singularity, but 
the mass term repels from it. 

While it is possible to check various isolated cases, treating all the in- 
finitely many theories in a perturbative manner is not viable. There is the 
further drawback, which takes a disjunctive form, depending on the choice 
of field variable. In the w-adapted variables 7, the gravitational field means 
different things for different theories. On the other hand, if one uses the 
neutral w = field, then all theories give an infinite series — no case is a 
polynomial, making the analysis difficult. 

Fortunately one can avoid perturbative treatments altogether and extrem- 
ize the mass-interaction part of the Lagrangian with respect to \J—g. One 
readily finds that the only critical point is the expected vacuum \J—g = a/— r) 
(which gives 7 = for all w) and that it is indeed the ground state for all w. 
In this sense massive scalar gravity is stable for all values of w. Some of the 
theories repel infinitely from the singularity yp^g = 0, while others do not. 



7 Why Canonical Tensor Derivation Is Sim- 
ple Only for w = | Theory 

It is not difficult to see why Freund and Nambu discovered in effect the w — ^ 
theory but not any of the other massive scalar gravities found here. They 
use the canonical stress-energy tensor in its standard simple form, as in their 
equation 3b, where the trace is given as 

8C • -4£. 



It is well known that one is permitted to add terms with automatically van- 
ishing divergence to the stress-energy tensor, terms sometimes called "curls" 
[38] by virtue of their resemblance to the vector calculus theorem that the 
divergence of a curl is (itself a consequence of 1 — 1 = 0). It is quite 
understandable that in a 3-page paper Freund and Nambu did not con- 
sider this option (though terms something like this were studied, still in 
a non-gravitational context, by Mack, by Chang and Freund, and by Aurilia 
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[551 EHJ ED])- The derivation of the above massive scalar gravities using the 
canonical tensor for arbitrary w, one can show, requires in the trace of the 
canonical stress-energy tensor the d'Alembertian of the term 

7 (1 + 8w\ZiiG^) 2^ 1 

as expressed in Cartesian coordinates, the use of which id advantageous when 
the canonical stress-energy tensor is employed. This extra term, which has 
second derivatives, vanishes if and only if w — |. Thus the neglect of this 
term will cause the derivation to fail except in the case w — |. 

The scalar gravity theory of Dehnen and Frommert [26] is equivalent to 
the Freund-Nambu theory |24j if one restricts the latter to matter fields with 
conformally invariant kinetic terms. This class includes not only electromag- 
netism and Yang-Mills theories (spin 1), but also fermions (spin |). It does 
not include standard scalar fields, however, though they do not remark on 
that important limitation. With (p being gravity and \ being matter, Fre- 
und and Nambu find that a standard scalar field coupled to gravity has an 
interaction involving <f)(dx) 2 - Dehnen and Frommert assume without justifi- 
cation that no such terms exist, as well as assuming that no nonlinear terms 
<p(d(f)) 2 arise. The latter terms can be absorbed by a nonlinear field redefi- 
nition of the gravitational potential <fi. For electromagnetism and Yang-Mills 
fields, taking the potentials to be covectors (as one usually does) removes 
the (p(dx) 2 terms. For spinors, I recall from above that redefining the spinor 
fields to have density weight | removes the coupling of the gravitational po- 
tential to the spinor (in this term roughly like (fnfjdi/j). The ability of 
the Dehnen-Frommert theory to accommodate spin | and spin 1 fields was 
exploited in some subsequent papers where only those matter fields were en- 
tertained j63j [TTJ . The Dehnen-Frommert derivation is not so simple because 
it involves multiple field redefinitions motivated by the need to recover from 
the assumption of purely nonderivative coupling between gravity and matter 
or the desire to derive massive scalar gravity in a Higgs-looking fashion. 

In view of the disadvantages of the (unsymmetrized, unimproved) canon- 
ical energy-momentum tensor for various fields — asymmetry, gauge depen- 
dence, and nonvanishing trace in some contexts where one might have wanted 
the trace to vanish [72] — the use of the metrical definition is a good deal more 
convenient. In principle one could use a Belinfante-Rosenfeld equivalence the- 
orem and employ the symmetric Belinfante tensor or the like. However, it 
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seems not very practical to derive an unknown Lagrangian density for non- 
scalar matter fields by looking for solutions to a messy identity involving all 
sorts of partial derivatives of the Lagrangian density with respect to the field 
derivatives. It seems not accidental that thus far only via the metrical defi- 
nition has an infinity, or even a variety, of universally coupled scalar gravities 
been derived. 

As Josep Pons has recently recalled [73], the process for making a flat 
space-time theory formally generally covariant admits considerable freedom 
in choosing density weights for fields: the comma-goes-to-semicolon rule thus 
has considerable ambiguity which tends to go unnoticed. While for many pur- 
poses the choice of density weight does not matter much (other than affecting 
the forms of the Lie and covariant derivatives) , the metrical stress-energy ten- 
sor, in particular its trace, is significantly affected. Above I observed that the 
use of a suitably densitized spinor allowed y/^g to disappear completely from 
the massless Dirac equation; I observe that one can do the same thing using 
conformally invariantly coupled scalar fields. The scalar field is replaced by 
a scalar density <p w of weight w = (which comes to \ in four space-time 
dimensions) by absorbing suitable powers of \/—g; the pleasant result is that 
\J—g thereupon disappears completely from the theory. The expression 




is the same for all conformally related metrics, from which it follows that 



\f—g~ simply cancels out altogether. Multiplying this expression by 4> w gives 

a scalar density of weight 1, which is thus a suitable Lagrangian density. 
One should be able to expand the expression out using g^ u and a/~# an d 
watch the latter disappear; the result will not look like a scalar density, so the 

2 

expression y/—g~[V 2 4> w — R<fiw) has some advantages. One could also 
drop some total divergences to remove second derivatives of (fi w and/or the 
metric, if desired, perhaps at the expense of manifest covariance (somewhat 
like the Einstein rr Lagrangian density for General Relativity). The absence 
of \J —g immediately implies a traceless metrical stress-energy tensor even off- 
shell. The use of densitized scalar and spinor fields thus allows one to identify 
and reject \J— g as surplus structure in some notable contexts. 
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